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The K-fi Shadowed Fading Model: 
Unifying the K-fi and rj-ji Distributions 

Laureano Moreno-Pozas, F. Javier Lopez-Martinez, Jose F. Paris and Eduardo Martos-Naya 


Abstract —This paper shows that the recently proposed k-/i 
shadowed fading model includes, besides the model, the 
Tj-fi fading model as a particular case. This has important 
relevance in practice, as it allows for the unification of these 
popular fading distributions through a more general, yet equally 
tractable, model. The convenience of new underlying physical 
models is discussed. Then, we derive simple and novel closed-form 
expressions for the asymptotic ergodic capacity in K-fi shadowed 
fading channels, which illustrate the effects of the different 
fading parameters on the system performance. By exploiting the 
unification here unveiled, the asymptotic capacity expressions for 
the K-fi and rj-fi fading models are also obtained in closed-form 
as special cases. 


I. Introduction 


The scientific community has been recently interested in the 
definition of new generalized fading models, aiming to provide 
a better fit to real measurements observed in different scenarios 
. In such context, the K-/i and r/-/r fading models have 
become very popular in the literature due to their versatility 
to accommodate to different propagation conditions and their 
relatively simple tractable form Ih MlOl] . 

The and 77-/1 fading models, first introduced in El 
and I 12 I] . were independently derived to characterize very 
different propagation conditions. On the one hand, the k- 
ji distribution can be regarded as a generalization of the 
classic Rician fading model for line-of-sight (LOS) scenar¬ 
ios, extensively used in spatially homogeneous propagation 
environments. On the other hand, the 77-/7 distribution can 
be considered as a generalization of the classic Nakagami- 
q (Hoyt) fading model for non-LOS scenarios, often used in 
non-homogeneous environments. Therefore, and because they 
arise from different underlying physical models, there is no 
clear connection between the K-/t and 77-/7 fading models. 

One of the most appealing properties of the k -/7 and 
77-/7 fading models is that they include most popular fading 
distributions as particular cases. For instance, the Rician, 
Nakagami- 777 , Rayleigh and one-sided Gaussian models can 
be derived from the k -/7 fading by setting the parameters k 
and /7 to specific real positive values. Similarly, the 77-/7 fading 
model includes the Nakagami-g, Nakagami-777, Rayleigh and 
one-sided Gaussian as special cases. 
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Very recently, the 7 c -/7 shadowed fading model was intro¬ 
duced in || 3 l, with the aim of jointly including large-scale 
and small-scale propagation effects. This new model exhibits 
excellent agreement when compared to measured land-mobile 
satellite EjIi . underwater acoustic II Gl and body commu¬ 
nications fading channels 1 15 , 113 ], by considering that the 
dominant components are affected by random fluctuations Et] . 
This model | ^includes the popular Rician shadowed fading 
distribution 11311 as a particular case, and obviously it also 
includes the 7 c -/7 fading distribution from which it originates. 
However, as we will later see, the versatility of the k -/7 
shadowed fading model has not been exploited to the full 
extent possible. 

In this paper we show that the k -/7 shadowed distribution 
unifies the set of homogeneous fading models associated with 
the K -/7 distribution, and strikingly, it also unifies the set 
of non-homogeneous fading models associated with the 77-/7 
distribution, which may seem counterintuitive at first glance. 
In addition to a formal mathematical proof of how the main 
probability functions introduced by Yacoub originate from the 
ones derived in Ell, we also establish new underlying physical 
models for the k -/7 shadowed distribution that justify these 
phenomena. In fact, we propose a novel method to derive the 
Nakagami-g (Hoyt) and the 77-/7 distributions which consists in 
using the shadowing of the dominant components to recreate a 
non-homogeneous propagation environment. This connection, 
which is here proposed for the first time in the literature, 
has important implications in practice: first, and contrary to 
the common belief, it shows that the k -/7 and 77-/7 fading 
distributions are connected. Hence we can jointly study the k -/7 
and 77-/7 fading models by using a common approach instead 
of separately E. 13 . iSj. Besides, it implies that when deriving 
any performance metric for the k -/7 shadowed fading model, 
we are actually solving the same problem for the simpler k -/7 
and 77-/7 distributions at no extra cost. 

Leveraging our novel approach, we derive simple and 
closed-form asymptotic expressions for the ergodic capacity of 
communication systems operating under k -/7 shadowed fading 
in the high signal-to-noise ratio regime, which can be evidently 
employed for the k -/7 and 77-/7 distributions. Unlike the exact 
analyses in El and S which require the use of the Meijer 
G- and bivariate Meijer G-functions, our results allow for a 
better insight into the effects of the fading parameters on the 
capacity. 

The remainder of this paper is structured as follows. In 
section II, we introduce the notation, as well as some defi¬ 
nitions and preliminary results. In section III, we propose new 
physical models for the k -/7 shadowed distribution. In section 
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IV we show how the K-fj, and rj-fj, distributions naturally 
arise as particular cases of the K-/r shadowed fading model. 
In section V, we use these results to investigate the ergodic 
capacity in shadowed fading channels and thus for the 
K-p, and p-p channels. In section VI, numerical results are 
presented. Finally conclusions are drawn. 


II. Notation and Preliminaries 


Throughout this paper, we differentiate the complex from 
the real random variables by adding them a tilde on top, so that 
X is a real random variable and 2 is a complex random variable. 
E[-] is the expectation operator. The symbol ^ signifies 
statistically distributed as while i.i.d means independent and 
identically distributed. Thus, x ^ Af{x, a^) symbolizes that 
X is statically distributed as a real Gaussian random variable 
with mean x and variance while z ^ CJ\f{z,a^) means 
that z is statistically distributed as a circularly symmetric 
complex Gaussian random variable with mean z and variance 
(T^. Moreover, although the classic fading distributions are 
usually defined in the literature by its envelope probability 
density function (pdf), we here present the results in terms of 
power probability density functions. 

Next, we present some definitions which are employed in 
the paper. 

Definition 1: Higher-order amount of fading M- 
The nth-order amount of fading is defined as 




E[7" 


- 1 


( 1 ) 


for the instantaneous SNR 7, where 7” = E[7]”. 

Definition 2: The generalized hypergeometric function. 

The generalized hypergeometric function of one scalar argu¬ 
ment is defined as 


+ 00 

pJ' q{(llf ■ - ■ , cip'-, ^ 1 , ■ - ■ ; bq] x) = ^ ^ 

r—0 


{ai)r ... {ap)r 

{bi)r . . . {hq)r r\ ’ 


where (a)r is the Pochhammer symbol eq. (6.1.22)], 
Oi G C and bj G C*\Z“. 

Definition 3: The gamma distribution. 

Let It; be a random variable which statistically follows a 
gamma distribution with shape parameter a and rate parameter 
P, i.e, w ^ r(Q;, P), then its pdf is given by 


fw{w) = 


P° 


r(a) 




(3) 


where r(-) is the gamma function 11221 eq. ( 6 . 1 . 1 )] and 

a,P G R+. 

Definition 4: The K-p, shadowed distribution j^]. 

Let 7 be a random variable which statistically follows a n-p 
shadowed distribution with mean 7 = E[7] and non-negative 


real shape parameters k, p and m, i.e, 7 ^ 
then its pdf is given by 




flil) = 


T{p)p{pK -\- m)"* 
X e iFi 




m\p] ■ 


+ k) 7 
pK m 7 


(4) 


where iFi{-) is the confluent hypergeometric function of 
scalar argument 0221 eq. (13.1.2)], which is a particular case 
of eq. (| 2 ]l. 

Definition 5: The K-p distribution fl. 

Let 7 be random variable which statistically follows a. K-p 
distribution with mean 7 = E[7] and non-negative real shape 


parameters k and p, i.e, 7 
by 


Ml) = 


p{l k) 


Sk.Yt'^ m)’ '^hen its pdf is given 


(5) 



where Ik{-) is the /c-th order modified Bessel function of 
first kind, which can be defined in terms of the Bessel 


hypergeometric function o-^i(’) 11221 eq. (9.6.47)1. 

Definition 6: The p-p distribution (Format 1) 111 . 

Let 7 be random variable which statistically follows an rj-p 
distribution with mean 7 = E[7] and non-negative real shape 


parameters 77 and p, i.e, 7 ■ 
by 


5^^ (7; 77 , p), then its pdf is given 


fih) = 


v^(l -f r]Y+ipf^+^ 

r(/i)7V^(i - 




X e 


2777 J 


/ p{l - pf 7 


277 


(6) 


III. Appropriate Physical Models 

In this section, we discuss the underlying physical model 
of the K-p shadowed presented in H. Then we propose two 
more general and suitable physical models which will allow 
us to later unify the K-p and 77 -p distributions. 


A. Previous model 

The K-p shadowed fading is presented as the generalization 
of the K-p model in H. Like in the K-p fading case, the signal 
can be decomposed into different clusters of waves, where 
each cluster has scattered waves with similar delays and the 
delay spreads of different clusters are relatively large J^]. In the 
K-p shadowed fading, the difference appears in the dominant 
component of each cluster, which is no longer deterministic 
and can randomly fluctuate because of shadowing. In fact, 
when we consider that the dominant component of each cluster 
suffers from the same shadowing, we can express the total 
signal power in terms of the in-phase and the quadrature 
components of each cluster as || 3 l 

w = + {y^ + ^q^Y, (7) 

i=l 

where p is the number of clusters; Xi and yi are independent 
real Gaussian random variables with zero mean and variance 
Y, which represent the real and imaginary parts of the i-th 
cluster scattering waves; pj and are the powers of the 
real and imaginary parts of the i-th dominant component, 
respectively, and ^ is the shadowing component which is 
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statistically distributed as a Nakagami-m random variable with 
shape parameter m and E[f^] = 1. 

It is known that this physical model follows a K-fi shadowed 
distribution ||3l, i-c., the instantaneous signal-to-noise ratio 
(SNR) 7 = jW/W, with 7 = E[7], W = E[kF] = d'^ + 
and = J2i=iPi + dh distributed as k, P,m), 

where the parameter k = /2a^ii represents the ratio be¬ 
tween the total power of the dominant components and the 
total power of the scattered waves. It is worth noticing that 
this result can be extended for /r taking non-integer positive 
values l^l, despite the model loses its physical meaning. 

However, this model imposes that the shadowing ^ is 
statistically distributed as a Nakagami-m random variable, 
which is a strict condition that is relaxed in the next section. 

B. Generalized model with the same shadowing for all the 
clusters 

The previous model in eq. dTJi clearly separates each cluster 
in the real and imaginary power components, so that the 
model can be defined by only using real random variables. 
Thus, using this time complex random variables, it can be 
reformulated as 

W = ( 8 ) 

i=l 

where 2 ) and 'pi can be related to the variables of the previous 
model in form of Zi = Xi -\- jpi and Pi = Pi + jqi- Hence, 
Zi ^ CAf{ 0 , 2cr^) represents the scattering wave of the i-th 
cluster and is the deterministic dominant power of the 
i-th cluster. 

A straightforward generalization of the previous model is to 
consider a complex shadowing component, so that we obtain 
the following model 

(9) 

i=l 

where ^ is now a complex random variable, with 
^ r(TO,m) and arbitrary phase. 

This new model obviously represents a similar scenario as 
the previous one in Section III.A, since all the clusters suffer 
from the same shadowing which can be justified by the fact 
that the shadowing can occur near the transmitter or receiver 
side. The pdf of the instantaneous SNR of the model in eq. dUl 
is derived as follows. 

Lemma 1: Let 71 = fiLli/Cli, with 71 = E[7i]> 

be the instantaneous SNR of the model in eq. (|9l). Then 
7i ^ where k = J2i=i IPiP/Scr^/r. 

Proof: The signal power Hi conditioned to the shadow¬ 
ing power follows a K-p distribution. Moreover, since 

then kP Eti ^ r(m, mEti 
Thus, the conditional form here obtained is equivalent to the 
one in jst] and so we can follow the same steps to prove that 
the SNR of the model follows a K-p shadowed distribution. 

In fact, we have proved that the distribution of the model 
is independent of the phase of the shadowing component 
In the next section we propose another physical model for the 
K-p shadowed distribution. 


C. Generalized jnodel with Ltd. shadowing 

When the shadowing does not occur near the transmitter 
or receiver sides, all the clusters could suffer from different 
shadowing effects, so the instantaneous received power can be 
expressed as 

( 10 ) 

i=l 

where 2) ^ CAf{0, 2cr^); |pp is the power of each dominant 
component; and is a complex random variable which 
represents the shadowing component of the i-th cluster, where 
are i.i.d., i.e, ^ r(m,m) Vi. 

Actually, these propagation conditions are likely to occur 
in real scenarios where the dominant components of different 
clusters could travel through different paths that are separated 
enough to suffer from independent large-scale propagation 
effects. The pdf of the instantaneous SNR of the model in 
eq. ([Tol l is derived next. 

Lemma 2: Let 72 = 72 H 2 /H 2 , with 72 = E[ 72 ], be 
the instantaneous SNR of the model in eq. (fTOl i. Then 
72 tn) where k = |pp/ 2 cr^/r and m = p-m. 

Proof: The conditioned signal power H 2 conditioned to 
the sum of the i.i.d. shadowed dominant component powers 
P — follows a K-p distribution. Moreover, 

since ^ T{rh^m) Vi, then P ^ r(m, |ppm), where 
m = m. Thus, again we have the same conditional form 
as in jsl] and so we can follow the same steps to prove that 
the SNR of the model follows a K-p shadowed distribution. Q 
Therefore, the SNRs of both physical models presented in 
eq. (|9]l and eq. (fTOl l follow a K-p shadowed distribution. The 
closed-form expressions for the cumulative distribution and 
the moment generating functions can be found in fl. 

IV. K-p AND p-p Unification 
In the previous section, we have introduced two different 
physical models which lead to the K-p shadowed distribution. 
Now, we show how each of these models reduces to the general 
K-p and p-p fading distributions, respectively. By doing so, we 
show that the K-p shadowed distribution can unify all classic 
fading models, both for homogeneous and non-homogeneous 
propagation conditions, and their most general counterparts 

A. K-p distribution and particular cases 

The K-p distribution is destined to model homogeneous 
environments, where the scattering for each cluster can be 
modeled with a circularly symmetric random variable. The 
derivation of the K-p distribution from Lemma 1 is given in 
the following corollary. 

Corollary 7: Let 71 = 71 H 1 /H 1 , with 71 = E[ 7 i], 
be the instantaneous SNR of the model in eq. i.e., 

71 ^ K, p, m). If m ^ 00 , Pi ^ 5„^(7i; p)- 

^The different shadowing components do not have to be identically 
distributed to complete the proof. All that is needed is that the normalized rate 
parameter of each shadowing power must be equal Vf. Although 

from a mathematical point of view this is a valid model, this scenario is hard 
to imagine in practical conditions. Therefore, we will restrict ourselves to the 
case with i.i.d. shadowing components. 


4 


Proof : By taking the limit m —oo in eq. (IH and applying 
the following properties 

lim iJ^i(a-b;-z) = oJ^i( b;z), (11) 

lim + =e~^, (12) 

a->oo Y a J 

where the eq. (fT^ is the well-known limit that defines the 
exponential function, we obtain the pdf in eq. (| 5 ]l. 

Corollary 1 is interpreted as follows; the K-/i distribution 
is derived by eliminating completely the shadowing of each 
dominant component, which can be done by taking to —00, 
so that the dominant component of each cluster becomes 
deterministic. Actually, as the parameter to grows, the pdf of 
each dominant component is gradually compressed and, at the 
limit TO —4 00, it becomes a Dirac delta function. Thus, the 
model is defined by a circularly symmetric complex random 
variable with some non-zero mean in each cluster, so that we 
obtain the n-fi model, whereas in case that p, = 1 we have the 
Rician fading model. 

Next, we derive the Nakagami-m physical model from the 
physical model in eq. 

Corollary 2 : Let 71 = fiCli/Cli, with 71 = E[7 i], 
be the instantaneous SNR of the model in eq. dUl, i.e., 
71 ^ p, m). If K 0 , 71 r(p, p). 

Proof: By taking the limit k —0 in eq. (| 4 ]i and applying 
the following property 

lim pJ^q{ai.. .ap\bi.. .bq-,cz) = 1 , ( 13 ) 

c^O 

we obtain the pdf in eq. OJ. Notice that eq. (El can be 
carried out by simply exploiting the series expression of the 
hypergeometric function of scalar argument, where the first 
term has the unit value and the rest of the terms are powers 
of the scalar argument E eq. ( 13 . 1 . 2 )], so that they become 
zero when taking the limit. 

We give the following interpretation about Corollary 2 . By 
tending k — 0, we eliminate all the dominant components of 
the model, regardless the value of the shadowing parameter to, 
so that we only have scattering components in each cluster, i.e., 
we obtain a model which follows a Nakagami-TO distribution 
or one of its particular cases, Rayleigh or one-sided Gaussian, 
depending on the value of p. 

B. rj-fj, distribution and particular cases 

The Nakagami-g (Hoyt) and the 77-p distributions are em¬ 
ployed in non-homogeneous propagation conditions environ¬ 
ments, where the scattering model is non-uniform and can be 
modeled by elliptical (or non-circularly symmetric) random 
variables. At hrst glance, such scenario does not seem to 
ht with the K-p shadowed fading model. However, we can 
give a different interpretation to the cluster components of 
the physical model in eq. (EOl l: they can be interpreted as a 
set of uniform scattering waves with random averages. These 
random fluctuations in the average, which are different for each 
cluster, are responsible for modeling the non-homogeneity of 
the environment and ultimately lead to breaking the circular 
symmetry of the scattering model. We must note that a 


similar connection was inferred in Il 23 h . where the Nakagami-q 
distribution was shown to behave as a Rayleigh distribution 
with randomly varying average power. We show next how the 
circular symmetry of the model can be broken by using the 
result of Lemma 2 . 

Corollary 3 : Let 72 = 72H2/H2, with 72 = E[72], 
be the instantaneous SNR of the model in eq. (EOl l. i.e. 


72 


5 , 


Kfimir, 1^1 p, to). If TO = p/ 2 , 72 


2 k +1 ’ A^/2)' 

Proof: When to = p/ 2 , we can apply in eq. (| 4 |l the 
following property lE 3 . eq. ( 9 . 6 . 47 )] 


(a; 2a; z) = ("a + M (|) , 

( 14 ) 

and so we have the eq. (01 after some simple algebraic 
manipulations. 

Hence, we have shown that the 77-p fading distribution arises 
as a particular case of the more general K-p shadowed model. 
This is one of the main results in this paper. Notice that when 
TO = p /2 = 0.5 we obtain the Nakagami-g model with shape 

parameter q = -sj 2^+1 7 = 7^ for 7‘M format 1 jillFI 

We can give the following interpretation of the result in 
Corollary 3 . Let us consider that each shadowing component 
can be expressed as = Xi ■ e^‘^, where Xi ^ A/'( 0 , 1 ), so 
that ^ r(l/ 2 , 1 / 2 }. When we consider the simplest case 
on which the phase of is deterministic and set to zero, then 
becomes a real Gaussian random variable. Thus, we are 
adding a real Gaussian random variable to a complex Gaussian 
random variable, so that the circular symmetry of the model 
is broken. In the general case of an arbitrary phase for 
the circular symmetry would be broken in a direction of the 
complex plane different to the real axis. Therefore, while the 
K-p model in Section IV.A is obtained by totally eliminating 
the randomness of the shadowing component, this is not the 
case for the 77-p fading model. 

In turn, the Nakagami-m model can be also deduced from 
the K-p shadowed physical model of eq. ( EOt with a similar 
method, i.e., without eliminating directly the dominant com¬ 
ponent. 

Corollary 4 : Let 72 = 72H2/H2, with 72 = E[72], 
be the instantaneous SNR of the model in eq. (EOl) . i.e, 
72 ^ 5 «/xm( 7 ; ft, m). If TO = p, 72 -- r(p, p). 

Proof : The result is straightforward by applying 

iTi{a]a]z'^ = ( 15 ) 


and making some algebraic manipulations. 

Notice that by setting to = p, we transform the i.i.d. 
random dominant components of the eq. (EOl) into scattering 
components. In fact, since to = p • tti, then we set 771 = 1 
and the i-th random dominant component becomes a Gaussian 
random variable. Thus we are adding two Gaussian random 
variables together in each cluster, which leads to an equivalent 
Gaussian random variable, so that the one-sided Gaussian, 
Rayleigh or Nakagami-m models are obtained depending on 
the number of clusters p considered. 


^The rj-fi fading model (format 1 ) is symmetrical for 77 S [ 0 , 1 ] and rj S 
[1, 00 ]. We have q = ^/rj 01 q = depending on the interval. 
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TABLE I 

Classic and generalized models derived from the K-fi 

SHADOWED FADING 


Channels 

Shadowed Parameters 

One-sided Gaussian 

a.) = 0.5, /s —^ 0 

b) /T = 0.5, m = 0.5 

Rayleigh 

a) /i = 1, K ^ 0 

b) /j, = 1, m = 1 

Nakagami-m 

a.) = m, IS 0 

b) = m, rn = m 

Nakagami-g (Hoyt) 

ft = 1, K = {1 — q^)j2q^, m = 0.5 

Rician with parameter K 

^ = 1, K = K, m ^ oo 

K,-fl 

fj, = {i, K = K,, rn^oo 


= 2/i, K = (1 — r))l2r}, m = fL 

Rician shadowed 

^ = 1, ^ = K, m = m 


The table I summarizes all the models that are derived from 
the K-/r shadowed fading model, where the K-fi shadowed 
model parameters are underlined for the sake of clarity. When 
the K-fi shadowed parameters are fixed to some specific 
real positive values or tend to some specific limits, we can 
obtain all the classic central models, i.e., the Rayleigh, one¬ 
sided Gaussian, Nakagami-g and Nakagami-m, the classic 
noncentral Rician fading, and their general counterparts, the 
Rician shadowed, K-fi and rf-fi fading models. 

It is remarkable that there are two ways for deriving the one¬ 
sided Gaussian, Rayleigh and Nakagami-m models, depending 
on whether the approaches in Section IV. A or Section IV.B are 
used. 


V. Application: Ergodic Capacity in the High-SNR 
Regime 

The characterization of the ergodic channel capacity in 
fading channels, defined as 


C [bps/Hz] = 


^+oo 


log 2 (l + 7 )/ 7 ( 7 )d 7 , 


( 16 ) 


where 7 is the instantaneous SNR at the receiver side, has been 
a matter of interest for many years ll 24 l 427 ll . While for the case 
of Rayleigh fading it is possible to obtain relatively simple 
closed-form expressions for the capacity, the consideration 
of more general fading models | 19 , 20 , | 23 l leads to very 
complicated expressions that usually require the use of Meijer 
G-functions. 

In order to overcome the limitation of the exact characteriza¬ 
tion of K-fi shadowed channel capacity due to its complicated 
closed-form 0, it seems more convenient to analyze the 
high-SNR regime. In this situation, the ergodic capacity can 
be approximated by 112 iL eq. (8)] 


which is asymptotically exact and where L is a constant value 
independent of the average SNR that can be given by 


L = 



( 18 ) 


In fact, the parameter L can be interpreted as the capacity 
loss with respect to the additive white Gaussian noise (AWGN) 
case, since the presence of fading causes L > 0 . When there is 
no fading, L = 0 and this reduces to the well-known Shannon 
result. Using this approach, we derive a simple closed-form 
expression for the asymptotic capacity of the K-fi shadowed 
model, which is a new result in the literature. 

Lemma 3 : In the high-SNR regime, the ergodic capacity of 
a K-fi shadowed channel can be accurately lower-bounded by 


(7)177 “ 1082(7) 


( 19 ) 


where log2(-) is the binary logarithm, e is the base of the 
natural logarithm, 7 is the average SNR at the receiver side, 
i.e. 7 = E[7], and LK/im can be expressed as 




fiK-\- m 


k), 


■log 2 (e) 


K{fi — m) 


fiK - 


( 20 ) 


X.7 T'. 


( 1 , 1 ,/r-m-f l;2,p-f 1; —^ 

V fiKmi 


where '(/;(•) is the digamma function 0221 eq. ( 6 . 3 . 1 )] and 
3J^2(-) is a generalized hypergeometric function of one scalar 
argument. 

Proof-. See Appendix A. 

Notice that when p, = 1 , we obtain the ergodic capacity of 
the Rician shadowed in the high-SNR re gim e. 

As opposed to the exact analysis in 1201 ] . which requires 
for the evaluation of a bivariate Meijer G-function, Lemma 3 
provides a very simple closed-form expression for the capacity 
in the high-SNR regime. More interestingly, since the K-fi and 
rf-fi fading channel models are but particular cases of the K-fi 
shadowed distribution, we also obtain the capacity in these 
scenarios without the need of evaluating Meijer G-functions 
as in id. This is formally stated in the following corollaries. 

Corollary 5 : In the high-SNR regime, the ergodic capacity 
of a K-fi channel can be accurately lower-bounded by 


(7)177 = 1082(7) - L, 




( 21 ) 


where can be expressed as 

= - log2(e)'(/'(p) + log2(p) + log2(l + k ) 

g X / \ ( 22 ) 

- Klog2(e)2J^2(^l,l;2,p+ 1 ;-pkJ. 

Proof-. The eq. (l 22 l i is derived by applying the limit m —00 
in eq. (l20l) . so that the 3J^2(’) collapses in a 2-^2(■) hypergeo¬ 
metric function since 

2; 

lim 3-^2(ar,a2,c;6i,62;-) = 2-E2(ar, 02; 61,62; z)- ( 23 ) 

c —>^00 c 

Corollary 6 : In the high-SNR regime, the ergodic capacity 
of an Tf-fi channel can be accurately lower-bounded by 

^»?A‘(7)l77 ~ 1082(7) ~ Lrifi, 


C(7)l77 = 1082(7) - L, 


( 17 ) 


( 24 ) 
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TABLE II 

Ergodic capacity loss in the high-SNR regime eor dieeerent 

CHANNELS 


Channels 

Ergodic capacity loss (L) [bps/Hz] 

One-sided Gaussian 

1 -|-7e ■ log 2 (e) Ri 1.83 

Rayleigh 

7 e ■ log 2 (e) Ri 0.83 

Nakagami-m 

log2(m) - log2(e)i/)(m) 

Nakagami-g (Hoyt) 

1 -|-7e ■ log 2 (e) -|-log 2 ( 

Rician with parameter K 

log 2 (l-ll/A')-log 2 (e)r( 0 ,i^) 


- log2(e)i/’(M) + log2(M) + log2(l + «) 
-Klog2(e)2J2 (l, 1; 2, -1- 1 ; —ajk) 


-log2(e)V’(2Ai) +log2(M) 

+ log 2 (1 -1- »?) -1- log 2 (e) 
X3J'2(l,l,/i-yl;2,2Ai-yl;l-77) 

Rician shadowed 

7e ■ log2(e) 

log2 ( ,„(1+K) ) + (e) li+m 

X3-T2(i, 1,2 m;2,2;^^^) 


where can be expressed as 

= - log2(e)i/i(2/r) + + log2(l + ry) 

+ log 2 (e)^^—^ 3 -^ 2 (l, 1, M + 1; 2, 2/r + 1; 1 - ry). 

(25) 


Proof : We obtain the ly-yr asympototic capacity loss from 
eq. (|20] | by setting m= yr, yt= 2yr and k.= 2^ as Table I 
indicates. 

Hence, the expressions of the n-ji and ?y-/r asymptotic ca¬ 
pacities have been jointly deduced from the result of Lemma 3, 
which are also new results. Moreover, deriving the asymptotic 
capacity of the it-p shadowed has not been harder than 
deriving the K-/i or ry-p asymptotic capacities directly, since 
the K-fj, and rj-fx moments are expressed, like in the K-yr 
shadowed case, in terms of a Gauss hypergeometric function 
llill . Thus, we are hitting two (actually three) birds with one 
stone. 

Using the equivalences in Table I, we can obtain even 
simpler expressions for classic fading models which reduce to 
existing results in the literature, for Nakagami-m a , Rician 
128] and Hoyt |23, For the sake of clarity, we omit the 


straightforward derivations of the rest of asymptotic capacities. 
Instead, we summarize in Table II their capacity losses with 
respect to the AWGN channel in the high-SNR regime, where 
r(a, b) is the incomplete gamma function ll22[ eq. (6.5.3)] and 
7 e is the Euler-Mascheroni constant, i.e., 7 e ~ 0.5772. 


VI. Numerical Results 

We now study the evolution of the capacity loss for the 
K-yt shadowed, K-yt and ry-yi fading models with respect to the 
AWGN case. In Fig.[T]and Fig.|2] we plot the ergodic capacity 
of the classic and generalized fading models, respectively. 


2 

'O 



Fig. 1. Comparison of classic channel ergodic capacities with their asymp¬ 
totic values in the high-SNR regime. 


5 

c/5 

Cu 


1^ 

'O 



7 [dB] 


Fig. 2. Comparison of generalized channel ergodic capacities with their 
asymptotic values in the high-SNR regime. 


We observe that all the models converge accurately to their 
asymptotic capacity values, remaining below the Shannon 
limit, i.e, the capacity of the AWGN channel. Therefore, the 
asymptotic ergodic capacity expression derived in Lemma 3 
for the K-/i shadowed model is here validated for the one¬ 
sided Gaussian, Rayleigh, Nakagami-m, Nakagami-g, Rician, 
Rician shadowed, K-yi and ly-yt ergodic capacities in the high- 
SNR regime. 

In Figs. [316] we show the evolution of the it-yi shadowed 
asymptotic capacity loss when m grows. When the shadowing 
cannot be negligible, i.e, for m < 3 which corresponds to 
Figs. [3|5] having more power in the dominant components 
does not always improve the ergodic capacity, but sometimes 
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Fig. 3. Evolution of the shado\ved ergodic capacity loss in the high-SNR 
regime for fixed m = 0.5. 


raises considerably the capacity loss, especially for a great 
number of clusters. When m > 20, the shadowing can be 
neglected and the model actually tends to the K-/r fading, 
where an increase in the power of the dominant components 
is obviously favorable for the channel capacity. Therefore, 
receiving more power through the dominant components does 
not always increase the capacity in the presence of shadowing. J- 
In fact, we observe two different behaviors in the capacity 
loss evolution with respect to the parameter k. For m < fjL, 
increasing the parameter k is detrimental for the capacity. 
Conversely, when to > /i the capacity is improved as k is 
increased, i.e. in the presence of a stronger LOS component. 

In the limit case of to = /r, we see that the capacity loss is 
independent of k. We also see that the capacity loss decreases 
as /r grows, since having a larger number of clusters reduces 
the fading severity of the small-scale propagation effects. 

We have also marked in Figs. [36] some models that can be 
deduced from the k-^ shadowed model. In fact, we can see 
them in the different legends and also at some specihc points 
rounded by a circle in different curves. 

Finally, Fig. |7| and Fig. [8] depict the asymptotic ergodic 
capacity loss for the n-fi and rj-fj, fading models respectively. 

We observe that Fig. [7] is quite similar to Fig. | 6 ] because, as 
mentioned before, the K-fj, shadowed model with to > 20 
can be approximated by the k-h fading model. In Fig. 8 , 
we see that, regardless the number of clusters /r, there is a 
minimum in the channel capacity loss at ry = 1 which divides 
in two symmetric parts the fading behavior as expected. It also 
noticeable that in Fig. [ 8 ] we have specified the limit cases for 
p —^ 0 and p —^ oo. In fact, when /i = 0.5, the 7 y-/i collapses 
into the one-sided Gaussian model for 77 = 0 or p —^ 00 , 
whereas for 77 = 1 it collapses into the Rayleigh model. In 
turn, when /r = 1, the rj-fi is reduced to the Rayleigh case for 
77 = 0 or 77 00 . This is shown in the figure by including 

also the Rayleigh and one-sided Gaussian capacity loss values 



Fig. 4. Evolution of the k,-{i shado\ved ergodic capacity loss in the high-SNR 
regime for fixed m = 1. 



K 


Fig. 5. Evolution of the n-ii shado\ved ergodic capacity loss in the high-SNR 
regime for fixed m = 3. 

with horizontal dotted and dashed lines respectively. 

VII. Conclusions 

We have proved that the K-fi shadowed model unifies 
the K-/r and 77-/7 fading distributions. By a novel physical 
interpretation of the shadowing in the dominant components, 
we have shown that the k-/ 7 shadowed model can also be 
employed in non-homogeneous environments, which gives 
the 7C-/7 shadowed distribution a stronger flexibility to model 
different propagation conditions than other alternatives, when 
operating in wireless environments. Thus, the k-/ 7 shadowed 
model unifies all the classic fading models, i.e., the one-sided 
Gaussian, Rayleigh, Nakagami-TO, Nakagami-q and Rician 
fading channels, and their generalized counterparts, the k-/7. 






















































2 


2 



Fig. 6. Evolution of the K-fi shadowed ergodic capacity loss in the high-SNR 
regime for fixed m = 20. 



Fig. 8. Evolution of the rf-fi ergodic capacity loss in the high-SNR regime. 
The Rayleigh and the one-sided Gaussian capacity loss particular cases are 
also included in horizontal dotted and dashed lines respectively. 



Fig. 7. Evolution of the k,-{i ergodic capacity loss in the high-SNR regime. 


Appendix A 

Ergodic capacity of the shadowed fading in 
THE high-SNR regime 


In the high-SNR regime, it is well-known that the ergodic 
capacity can be lower-bounded by the eq. (HU), where the 
parameter is related with the nth-order derivative of 

the amount of fading , such as 1211 


= -log2ie)—AF!f 


n—0 


(26) 


Thus, thanks to eq. dB, we first have to compute the moments 
of the SNR at the receiver side before deriving EK/xm, i e. 


p-l-oo 

E[ 7 "] = / 

Jo 


+ mY' 

/• + 0O 






(27) 


Xi m,p; 


fYK{l + k ) 7 


fiK + m 7 


- dj. 


p-p, and Rician shadowed fading models. Using this connec¬ 
tion, simple new closed-form expressions have been deduced 
to evaluate the ergodic capacity in the high-SNR regime for the 
K-fj, shadowed, and hence for the n-p and p-p fading models, 
giving us clear insights into the contribution of the fading 
parameters on the capacity improvement or degradation. 


Observing that the remaining integral correspond to a Laplace 


transform evaluated in s = 
eq. (4.23.17)] 


At(l-tK) 


then have 


E[7"] = 


r(p + n) 7”m™(l -I- k) 

r(p) Y^{fjLK + mY 

X 2 F 1 ( m;p -f n;p; 


IJ.K 


F m 


(28) 


As a closing remark, one can think of whether the name of 
K-p shadowed distribution is still appropriate for this model, 
since its flexibility transcends the original characteristics pre¬ 
sented in i^j. 


where is the Gauss hypergeometric function of scalar 

argument 11221 eq. (15.1.1)]. 

By making a well-known transformation involving the argu¬ 
ments of the Gauss hypergeometric function 0221 eq. (15.3.3)], 
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we obtain 


r(/i) \^lm{l + K)) 

X 2^1 [ ^^■ — m, —n; ■ 




(29) 


HK + m^ 

The amount of fading is then deduced by using the product rule 
in eq. (|29] | with the Gauss hypergeometric function expressed 
in series form 11221 eq. (15.1.1)]. As the derivative of a 
Pochhammer symbol can be given by a difference of digamma 
functions we obtain 

d .pin) r(/x + n) / ^J,K + m Y 

i ^ rfw) Vumfl + K)/ 


dn 


r(/x) V + k) 

X I ?/>(/x + n) + log 
X 2-7^1 ~ tTj,—n;/i; 


fiK + m 
p,m(l + k) 

pK 

+ m 


+ 00 

-E 

r—l 


{fj. - m)r{-n)r 


\ / 

r\ 


i^^)r 

■y 


(y(—n + r) — jp{—n)) 


(30) 


where the infinite sum starts at r = 1 because the first term 
equals zero. Setting the moments order n = 0, we get 


n—0 


L,f.n,=-\og^{e)—AF(^^^ 

= - log,(e)V>(rt - log. ( y'jyy ) 

, / . ~ w)r(l)r-l inK+m) 

+ log2(e) 2^- 


( 31 ) 


(/x)r 


By applying some algebraic manipulations, we finally obtain 

= - log.(e)*(rt - log. ( 2)22 ) 

«:(/i — m) 


log^e)- 


fin ■ 


+ 00 

E 


(p - m + l)r i(l)r-l(l)r-l i^K+m) 


(P- T 1)1-—l(2)r—1 


(r- 1)! 


(32) 


where the infinite sum can be expressed in terms of the 
generalized hypergeometric function 3 Jy(-) and so we have 
the result of eq. (l20l) . Notice that this result gives a simple 
new expression for the derivative of the Gauss hypergeometric 
funcion c, z) with respect to a or b, when this 

same parameter a or 6 equals zero. In fact, this derivative is 
expressed in terms of the generalized hypergeometric function 
3 J^ 2 (’) instead of in terms of a Kampe de Feriet function as 
proposed in yj. |. 
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